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This  paper  and  Paper  I  in  this  series  (P.H.  Berens  and  K.R.  Wilson,  J.  Chem. 
Phys.,  in  press)  indicate  that  infrared  and  Raman  rotational  and  fundamental 
vibrational-rotational  spectra  of  dense  systems  (high  pressure  gases,  li¬ 
quids  and  solids)  are  essentially  classical,  in  that  they  can  be  computed 
and  understood  from  a  basically  classical  mechanical  viewpoint,  with  some 
caveats  for  features  in  which  anharmonicity  is  important,  such  as  the  de¬ 
tailed  shape  of  Q  branches.  It  is  demonstrated  here,  using  the  diatomic 
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20.  ABSTRACT  (continued) 

case  u  an  example,  that  ordinary,  i.e.  nonresonant,  Raman  band  contour*  can 
fce  computed  froa  classical  mechanics  plus  simple  quantua  corrections.  Classi¬ 
cal  versions  of  molecular  dynamics,  linear  response  theory  and  ensemble  aver¬ 
aging,  followed  by  straightforward  quantiss  corrections,  are  used  to-  compute 
the  pure  rotational  and  fundaaental  vibrational-rotational  Ranan  band  contours 
of  for  the  gas  phase  and  for  solutions  of  In  different  densities  of  gas 

phase  Ar  and  in  liquid  Ar.  The  evolution  is  seen  froa  multiple  peaked  line 
shapes  characteristic  of  free  rotation  in  the  gas  phase  to  single  peaks  charac 
teristic  of  hindered  rotation  in  the  liquid  phase.  Comparison  is  nade  with 
quantum  and  correspondence  principle  classical  gas  phase  spectral  calculations 
and  with  experimental  measurements  for  pure  Nj  and  in  liquid  Ar.  Three  ad¬ 
vantages  are  pointed  out  for  a  classical  approach  to  infrared  and  Raman  spme- 
tra.  First,  a  classical  approach  can  be  used  to  compute  the  spectra  of  com¬ 
plex  molecular  systems,  for  example  of  large  molecules,  clusters,  liquids, 
solutions,  and  solids.  Second,  this  classical  approach  can  be  extended  to 
compute  the  spectra  of  non-equilibrium  and  time-dependent  systems,  for  example 
infrared  and  Raman  spectra  during  the  course  of  chemical  reactions.  Third,  m 
classical  viewpoint  allows  experimental  infrared  and  Raman  spectra  to  be  un¬ 
derstood  and  interpreted  in  terms  of  atomic  motions  with  the  considerable  aid 
of  classical  models  and  of  our  well-developed  classical  intuition.  i 
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ABSTRACT 

This  paper  •an4-Paper  +  in  this  series  (P.  H.  Berens  and  K.  R.  Wilson,  J 
Lhem.  Pliys .,  in  press)  indicate  that  infrared  and  Raman  rotational  and  funda¬ 
mental  vibrational-rotational  spectra  of  dense  systems  (high  pressure  gases, 
liquids  and  solids)  are  essentially  classical,  in  that  they  can  be  computed  and 
understood  from  a  basically  classical  mechanical  viewpoint,  with  some  caveats 
for  features  in  which  anharmonicity  is  important,  such  as  the  detailed  shape  of 
0  branches.  It  is  demonstrated  here,  using  the  diatomic  case  as  an  example, 
that  ordinary,  i.e.  nonresonant,  Raman  band  contours  can  be  computed  from 
classical  mechanics  plus  simple  quantum  corrections.  Classical  versions  of 
molecular  dynamics,  linear  response  theory  and  ensemble  averaging,  followed 
by  straightforward  quantum  corrections,  are  used  to  compute  the  pure  rota¬ 
tional  and  fundamental  vibrational-rotational  Raman  band  contours  of  N\  for 
the  gas  phase  and  for  solutions  of  Ns  in  different  densities  of  gas  phase  Ar  and 
in  liquid  Ar.  The  evolution  is  seen  from  multiple  peaked  line  shapes  charac¬ 
teristic  of  free  rotation  in  the  gas  phase  to  single  peaks  characteristic  of  hin¬ 
dered  rotation  in  the  liquid  phase.  Comparison  is  made  with  quantum  and 
correspondence  principle  classical  gas  phase  spectral  calculations  and  with 
experimental  measurements  for  pure  N>  and  Ns  in  liquid  Ar  Three  advantages 
are  pointed  out  for  a  classical  approach  to  infrared  and  Raman  spectra  ./First,  a 
classical  approach  can  be  used  to  compute  the  spectra  of  complex  molecular 
systems,  for  example  of  large  molecules,  clusters,  liquids,  solutions,  and  solids 
Second,  this  classical  approach  can  be  extended  to  compute  the  spectra  of  non¬ 
equilibrium  and  lime-dependent  systems,  for  example  infrared  and  Raman 
spectra  during  the  course  of  chemical  reactions.  Third,  a  classical  viewpoint 
allows  experimental  infrared  and  Raman  spectra  to  be  understood  and  inter¬ 
preted  in  terms  of  atomic  motions  with  the  considerable  aid  ol  classical  models 
and  of  our  well-developed  classical  intuition 
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I.  INTRODUCTION 

The  properties  of  matter  and  the  processes  it  undergoes  are  determined  not  just  by  its 
average  structure  but  also  by  its  dynamics,  its  time  evolution.  Thus  to  understand  these  pro¬ 
perties  and  processes  we  would  like  a  window  through  which  to  observe  atomic  motions.  As 
Gordon  has  emphasized,1  infrared  and  Raman  spectra  intimately  reflect  both  intra-  and  inter- 
molecular  motions,  and  thus  in  principle  can  provide  such  a  window.  In  practice,  if  we  use 
quantum  mechanics  as  our  link  between  time  evolution  and  spectra,  we  are  limited  to  the 
analysis  of  rather  simple  cases  by  computational  complexity  and  by  our  inability  to  easily  con¬ 
ceptualize  in  a  quantum  framework.  Semiclassical  approaches  can  carry  us  further,  but  not  as 
yet  far  enough  to  handle  more  than  simple  systems.  Thus  a  generalized  classical  approach 
which  could  successfully  link  molecular  dynamics  and  infrared  and  Raman  spectra  would  be 
very  useful,  for  example  in  discovering  the  dynamics  underlying  liquid  properties,  or  in  follow¬ 
ing  chemical  change  in  condensed  phases. 

In  the  first  paper  ’  of  (his  series,  which  wc  will  henceforth  call  Paper  I,  we  demonstrated 
for  a  diatomic  molecule  that  pure  rotational  and  fundamental  vibrational-rotational  infrared 
band  contours  can  be  computed  from  classical  mechanics  plus  simple  quantum  corrections.  In 
this  paper  we  apply  the  same  classical  approach  of  molecular  dynamics,  linear  response,  and 
ensemble  averaging,  followed  by  simple  quantum  corrections,  to  compute  diatomic  pure  rota¬ 
tional  and  vibrational-rotational  Raman  spectra.  While  one  could  introduce  aspects  of  quantum 
mechanics  directly  into  the  dynamics,  and  thus  adopt  a  semiclassical  molecular  dynamics 
approach,1  1  12  we  illustrate  that  the  spectra  arc  so  very  nearly  classical  that  a  much  simpler 
approach  is  usually  quite  adequate,  one  in  which  the  introduction  of  quantum  mechanics  is 
delayed  until  the  end  and  then  the  final  classical  spectra  arc  quantum  corrected.  We  demon¬ 
strate  that  this  technique,  which  we  call  a  "Newtonian"  approach  to  emphasize  its  basically  clas¬ 
sical  nature,  allows  the  computation  of  diatomic  pure  rotational  and  fundamental  vibrational- 
rotational  Raman  spectral  band  contours  In  order  that  rigorous  tests  may  be  made  versus 
other  theoretical  approaches  and  experiment,  we  have  chosen  to  treat  a  molecular  system,  N;  in 
the  gas  phase  and  its  solution  in  Ar,  which  is  quite  simple  Nevertheless,  one  of  the  virtues  of 
the  Newtonian  approach  to  infrared  and  Raman  spectra  is  the  relative  ease  with  which  it  can  be 
extended  to  more  complex  systems,  such  as  molecular  clusters,  large  molecules,  liquids,  solu¬ 
tions,  and  solids,  as  well  as  to  non-equilibrium  and  time-dependent  systems11  In  order  to 
demonstrate  (his  approach  and  to  develop  the  necessary  quantum  corrections,  we  compare 
theoretical  spectra  computed  in  three  different  ways  /)  quantum  mechanics  without  linear 
response  theory,  ii)  correspondence  principle  classical  mechanics  without  linear  response 
theory,  te  the  limit  of  quantum  mechanics  as  Planck's  constant  approaches  zero,  and  mi  ordi¬ 
nary  Newtonian  classical  mechanics  with  classical  linear  response  theory  and  ensemble  averag¬ 
ing  from  classical  statistical  mechanics.  In  addition,  wc  compare  with  several  experimental 
spectra  for  pure  N:  and  for  N,  in  Ar 


i 
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II.  THEORETICAL  TECHNIQUES 

We  use  the  same  basic  Newtonian  theoretical  approach  presented  in  Paper  I.  First,  classi¬ 
cal  mechanics  is  used  to  compute  the  time  evolution,  i.e.  the  atomic  trajectories,  of  the  system 
of  interest  from  a  set  of  initial  atomic  positions  and  velocities  and  a  given  potential  surface. 
Second,  the  time  varying  polarizability  matrix  for  the  entire  system  of  molecules  is  calculated 
from  the  atomic  trajectories.  Third,  linear  response  theory,  as  applied  to  Raman  spectra  by 
Gordon,1- 14  17  is  used  to  calculate  a  "specific”  Raman  spectrum  for  each  molecular  dynamics 
run.  Fourth,  these  "specific"  spectra  are  averaged  over  an  ensemble  appropriate  to  the  experi¬ 
mental  conditions  of  interest  through  the  choice  of  different  sets  of  initial  positions  and  veloci¬ 
ties  for  a  series  of  molecular  dynamics  runs.  Fifth,  quantum  corrections  are  applied  where 
necessary  to  the  purely  classical  spectrum  to  converge  toward  the  real  quantum  spectrum. 

In  Paper  I  we  showed  that  a  similar  set  of  five  steps  produce  surprisingly  accurate 
infrared  spectra.  We  show  here  that  Raman  spectra  are  also  well  reproduced  in  this  way,  using 
as  inputs  only  the  potential  energy  F(r)  and  the  polarizability  P(r)  as  functions  of  nuclear  posi¬ 
tions.  Since  the  main  aim  of  this  paper  is  to  show  that  relatively  simple  classical  techniques  and 
quantum  corrections  can  reproduce  quan  i  reality  for  Raman  spectra,  we  have  relegated  to  a 
series  of  appendices  most  of  the  complexity  of  the  quantum  treatment,  which  we  need  to  derive 
quantum  corrections  and  to  demonstrate  by  comparing  classical  and  quantum  spectra  that  a 
simpler  classical  approach  indeed  works.  These  append>es  can  be  skipped  by  those  readers 
only  interested  in  the  classical  results. 

A.  Molecular  dynamics  and  classical  trajectories 

Given  appropriate  inter-  and  intra-molecular  potentials  to  describe  the  interactions  among 
atoms,  molecular  dynamics  can  be  used  to  calculate  from  a  set  of  initial  positions  and  velocities 
the  time  evolution  of  the  system  of  interest  Computation  of  the  classical  trajectories 

r , ( z > . rN(r)  for  a  system  of  N  atoms  involves  numerically  integrating  the  -V  coupled 

differential  equations  given  by  Newton’s  Second  Law 

cl  F  r  ^‘r<  i  s'  mi 

--  - —  =-  F,  =  m.  — —  .  /  =  1 . A  ( 1 ) 

dr,  dr 

in  which  f  =  F(rt . rO  is  the  potential  energy  of  the  atoms  (or  more  accurately  the 

nuclei)  at  positions  r( . rv.  F,  =  F. (r| . rN)  is  the  force  on  the  /  th  atom,  and  m  is  its 

mass  A  modified  Verlet  integration  algorithm  is  used. Is- 14  For  the  case  below  of  an  N:  solu¬ 
tion  in  Ar,  minimum  image  periodic  boundary  conditions  with  truncated  octahedral  boun¬ 
daries’0  and  forces  smoothly  feathered  to  zero  over  the  outer  0.01  nm  of  the  radius  of  the 
inscribed  sphere  (which,  for  the  cases  considered  here,  is  always  greater  than  0.57  nm)  are  used 
to  reduce  edge  effects. 

This  type  of  calculation,  while  straightforward  in  concept,  can  be  extremely  demanding 
arithmetically.  We  have  built  tin  "instrument  for  theory",  partially  described  elsewhere.-'  which 
consists  of  a  generalized  program  package  and  a  network  of  computers  including  an  array  pro¬ 
cessor  and  a  three  dimensional  display  processor.  This  instrument  allows  us  to  compute  spec 
tral  and  other  properties  of  even  rather  large  and  complex  molecular  systems 

R.  Time  history  of  the  polarizability 

Given  the  Cartesian  tensor  P ( r , . r%)  which  describes  the  polarizability  in  space-fixed 

coordinates  of  the  system  as  a  function  of  atomic  positions,  we  use  the  trajectories  of  the 
atoms,  T|(r) . rN(r).  calculated  from  molecular  dynamics  to  compute  P ( / ) .  the  polarizabil¬ 

ity  tensor  of  the  system  as  a  function  of  time  Polarizability  tensors  in  molecule  or  bond  fixed 
coordinate  systems  are  transformed  into  space  fixed  coordinates  h>  the  use  of  direction 
cosines 


C.  Linear  response  and  specific  spectra 

Gordon11416  has  shown  that  linear  response  theory  (see  Appendices  A  and  B)  can  be 
used  to  compute  the  Raman  spectrum  by  relating  the  spectrum  o!'  fluctuations  which  the  polari¬ 
zability  tensor  P(r)  naturally  undergoes  whenRhe  system  is  free  of  outside  perturbation  to  the 
spectrum  of  the  Raman  scattered  light  which  Pit)  produces  when  the  system  is  driven  by  the 
oscillating  electric  field  of  light. 

From  Appendices  A  and  B,  the  appropriate  linear  response  equations  are 


dur 

dcudfl 

d2ir 

dindil 


7  <Tr(Pm,(0)  P„„(r)]> 


<Tr[P, 


in  which  the  left  hand  sides  of  Eqs.  (2)  and  (3)  are  the  differential  cross  sections  for  scattering 
into  angular  frequency  range  dcj  and  solid  angle  range  dCl  for  the  isotropic  and  anisotropic 
Raman  spectra, Respectively,  weighted  by  X,4  in  which  2ttX,  is  the  wavelength  of  the  scattered 
radiation.  The  P,„,  and  PlW,v,  are  the  isotropic  and  anisotropic  parts  of  the  polarizability  tensor 
of  the  system,  respectively,  as  defined  in  Eqs.  (A2)-(A4)  of  Appendix  A.  Tr  indicates  the  trace 
and  <  >  indicates  an  ensemble  average. 

As  Parseval’s  theorem33  and  the  Wiener-Khintchine  theorem17- 34- 3'  show,  the  Raman 
cross  sections  given  above  as  Fourier  transforms  of  time  correlation  functions  may  also  be  com¬ 
puted  in  the  mathematically  equivalent  form  of  power  spectra. 
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allowing  the  use  of  fast  Fourier  techniques.  Considerable  care  must  be  used  to  properly  apply 
spectral  estimation  theory,  windowing  and  windowing  corrections,3- 33- 36  in  both  the  correlation 
and  power  spectrum  methods  to  avoid  distorting  the  spectra  as  a  result  of  the  use  of  finite  time 
histories  Particular  care  must  be  taken  if  band  wings  are  to  be  correct  or  if  small  bands  arc  to 
be  seen  in  the  presence  of  large  ones.  We  use  a  four  term  —74  db  Blackman-Harris  window.36 
which  allows  the  correction  to  be  applied  as  a  simple  frequency-space  convolution  of  the 
Fourier  transform  of  the  polarizability  matrix  element  time  histories  with  the  Fourier  transform 
of  the  window  function,  in  a  parallel  manner  to  that  described  in  Paper  I  The  time  correlation 
functions  can  then  be  computed,  if  desired,  from  the  Fourier  inverses  of  Eqs.  <21  and  (3).  In 
this  way  spectral  estimation  theory  is  correctly  used,  undesirable  distortions  in  the  correlation 
functions  from  finite  time  histories  are  avoided,  and  the  speed  advantage  of  fast  Fourier  tech¬ 
niques  is  retained.34- 3> 

The  isotropic  and  anisotropic  cross  sections  may  be  related  to  the  usual  experimental 
parallel  and  perpendicular  spectra  by16  17  33 
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D.  Ensemble  averaging 

Since  a  specific  spectrum  is  computed  from  a  single  molecular  dynamics  run  representing 
only  a  small  portion  of  phase  space,  in  order  to  compute  a  spectrum  comparable  to  an  experi¬ 
mental  spectrum  taken,  for  example,  at  a  particular  temperature  and  pressure,  it  is  necessary  to 
average  the  specific  spectra  over  an  ensemble  of  runs  with  properly  chosen  different  initial  posi¬ 
tions  and  velocities.  To  provide  a  sufficiently  accurate  ensemble  average  we  must  insure  that 
enough  of  phase  space  is  sampled.  Thus,  to  begin  a  spectral  computation  we  choose  reasonable 

guesses  for  the  initial  positions  T|(0) . r %  (0)  and  equilibrate  the  system  to  the  temperature 

of  interest.  We  integrate  forward  in  time,  stopping  at  intervals  to  choose  a  new  set  of  velocities 
selected  randomly  from  a  Maxweil-Boltzmann  distribution.  After  the  system  is  sufficiently 
equilibrated,  we  begin  taking  specific  spectra.  The  final  positions  of  the  molecular  dynamics 
run  for  the  previous  specific  spectrum  are  used  as  the  initial  positions  for  the  next  run.  in  an 
attempt  to  sample  more  of  configuration  space  by  allowing  slow  rearrangements  in  configuration 
space  to  accumulate.  Because  for  weakly  interacting  systems  such  as  dilute  gases,  constants  of 
the  motion  for  individual  molecules  such  as  angular  momentum  are  conserved  over  long 
periods  (nonergodicity),  and  thus  even  long  molecular  dynamics  runs  will  sample  only  res¬ 
tricted  regions  of  phase  space,  it  is  necessary  to  introduce  some  stochastic  element  into  the 
deterministic  molecular  dynamics.  Therefore,  we  randomly  pick  a  new  set  of  velocities  from 
the  Maxweil-Boltzmann  distribution  at  the  beginning  of  each  run.  While  this  is  necessary  for 
weakly  interacting  systems,  it  may  contribute  little  to  strongly  interacting  ones  in  which  the 
velocity  autocorrelation  function  dies  off  quickly,  and,  by  turning  concerted  motion  through 
configuration  space  into  diffusive  motion,  it  might  even  slow  down  the  convergence  to  a  proper 
configuration  space  average  if  applied  too  frequently.14 

E.  Potential  and  polarizability  functions 

The  intramolecular  potential  function  used  for  N,  is  a  fourth  order  Taylor's  series  expan¬ 
sion  about  the  equilibrium  bond  distance-7  of  0.109768  nm.  The  potenna)  derivatives,  T",  I  . 
and  V'"\  at  the  equilibrium  bond  distance,  listed  in  Table  I,  are  derived  using  Eqs  (C21MC26) 
from  the  ground  stale  spectral  constants  listed  by  Huber  and  llerzberg.-1 


TABLE  I  Potential  function  parameters  for  N. 
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We  have  used  as  the  polarizabilities  for_N;  the  values  shown  in  Table  II.  where,  for  a  diatomic 
in  molecule  fixed  coordinates.  P  and  P,, are  the  parallel  and  perpendicular  components, 
respectively,  and  or,,,,  and  are  the  isotropic  and  anisotropic  polarizability  components, 

respectively,  given  by  Eqs  (A9)  and  ( A 1  () )  of  Appendix  A  The  linear  expansion  given  by 
Eqs.  (A29)  and  ( A 30)  is  used. 

TABLE  II.  Polarizabilities  for  N, 
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The  values  of  P" ,  P", ,  and  are  the  experimental  values  adopted  by  Morrison  and 
Hay  n  The  values  of  P\  ,  P‘,,.  and  o are  derived  from  the  value  for  after 

Asawarocngchai  ct  til.  *4  combined  with  the  value  of  the  absolute  differential  Raman  scattering 
cross  section  for  the  fundamental  vibrational-rotational  Q  branch  of  V  given  by  Schrottei  and 
Klockner1'  using  Eqs  ((16)  and  (Cl 7) 


F.  Quantum  corrections 

In  order  to  discover  what  corrections  are  necessary  to  bring  our  classical  calculations  into 
agreement  with  quantum  reality,  we  compare  in  Appendix  E  two  different  ways  of  computing 
the  Raman  spectrum:  /)  quantum  mechanics  without  linear  response  theory,  and  //)  correspon¬ 
dence  principle  classical  mechanics  without  linear  response  theory  as  given  by  the  limit  of  quan¬ 
tum  mechanics  as  Planck’s  constant  approaches  zero.  Our  Newtonian  classical  mechanics  with 
classical  linear  response  theory  and  ensemble  averaging  from  classical  statistical  mechanics 
should  agree  with  correspondence  principle  classical  mechanics  if  both  are  carried  out  accu¬ 
rately,  and  thus  if  we  learn  to  properly  quantum  correct  one  we  can  also  correct  the  other 

I.  Pure  rotational  bands 

The  detailed  balance  correction  to  the  pure  rotational  band  is  given  in  Eq.  <E5)  of 
Appendix  E  as 


.4  d'(T 
l'  dwd li 


;  /  h 
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dojdn 


=  exp(£7ta>/2) 


in  which  the  superscript  R  indicates  rotational  and  the  subscripts  Q  and  C  indicate  quantum  and 
classical,  respectively.  In  Eq.  (8),  j3=(Afl7T  1  in  which  kH  is  Boltzmann's  constant  and  T  the 
temperature.  This  correction  can  be  understood  by  considering  that  in  the  correspondence  prin¬ 
ciple  limit  as  7i—0  and  as  quantum  mechanics  converges  toward  classical  mechanics,  Raman 
scattering  at  frequency  w  corresponds  to  infinitesimally  small  quanta  of  energy  Ttu>  and  produces 
transitions  between  states  infinitesimally  close  in  energy,  and  thus  classically  of  essentially  the 
same  statistical  probability.  Quantum  mechanically,  lower  and  upper  slates  are  related  in  proba¬ 
bility  by  the  Boltzmann  factor  cxp(pltto),  and  thus  Stokes  and  Antistokes  cross  sections  must 
be  related  by  this  factor.  The  factor  of  one  half  arises  from  the  2tu  separating  the  Stokes  and 
the  Antistokes  angular  frequencies.  Thus  we  quantum  correct  a  pure  rotational  band  by  multi 
plying  it  bv  exp(£7ftu/2) . 

2.  t  ibrational-rotational  bands 

In  Appendix  D  we  make  the  approximation  of  separating  vibration  from  rotation,  and 
thus  we  compute  the  quantum  corrections  separately,  even  though  in  our  classical  molecular 
dynamics  we  make  no  such  separation.  The  quantum  correction  for  pure  harmonic  vibration,  as 
given  in  Eq.  (EIG)  of  Appendix  E,  is  an  intensity  correction 
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This  can  be  understood  by  considering  that  the  intensity  of  Raman  scattering  is  determined  by 
the  square  of  the  change  in  polarizability,  which  is  proportional  in  the  linear  approximation 
used  here  to  the  square  in  the  change  in  bond  distance,  and  thus  to  the  vibrational  energy 
Quantum  mechanically  the  oscillator  energy  for  tiuj»kHT  is  —If u>/2  while  classically  it  is 
much  lower.  — k,tT ,  and  thus  the  quantum  scattering  is  greater  than  the  classical  scattering 

It  is  also  shown  in  Appendix  E  that  a  relatively  small  frequency  shift  is  needed  in  order  to 
correct  for  anharmonicity  in  the  vibrational  potential  because  a  classical  oscillator  samples  the 
potential  near  the  bottom  of  the  well  at  —kHT  where  the  anharmonicity  is  smaller  than  in  the 
region  sampled  by  the  zero  point  and  higher  quantum  states.  Eor  the  cases  we  treat  here,  the 
anharmonic  corrections,  calculated  as  discussed  in  Appendix  E.  are  — 27.2  cm  1  at  83  K. 
-27.0  cm  1  at  96  K,  and  -23  6  cm  1  at  298  K.  The  cruder  approximation  -4wi',..v,.,  using 
spectral  constants  from  Huber  and  llcr/bcrg,’1  gives  -28  6  cm  '.  In  addition,  if  we  assume 
that  vibration  and  rotation  are  approximately  mechanically  and  thus  statistically  separable,  then 
the  same  detailed  balance  correction  given  by  Eq  (8)  needs  to  be  applied  within  each  vibra¬ 
tional  band,  this  time  multiplying  the  band  by  expi/iT?  (*}<«)/ 21.  in  which  Sm  is  measured  front 
the  rotationless  band  center 
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Therefore  the  quantum  correction  of  a  fundamental  vibrational-rotational  band  consists  of 
three  steps.  First,  the  anharmonic  ofTset  is  computed  by  the  method  given  in  Appendix  h  and 
the  band  is  shifted  to  its  quantum  corrected  position.  Second,  the  band  is  harmonically 
corrected  by  multiplying  by  the  factor  given  in  Fq.  (9).  Third,  the  vibrational-rotational  band  is 
corrected  for  detailed  balance  by  multiplication  by  exp[07?(Acu)/2l,  in  which  A u>  is  the  angular 
frequency  displacement  measured  from  the  rotationless  center  of  the  band 

III.  SPECTRA 

Fig.  1  shows  the  overall  band  contour  of  the  N:  dilute  gas  phase  Raman  spectrum  which 
we  will  now  examine  in  more  detail. 

A.  Gas  phase  rotational  spectra 

f  ig.  2  shows  the  anisotropic  pure  rotational  Raman  spectrum  for  gas  phase  NY  The  isotro¬ 
pic  rotational  Raman  spectrum  is  zero  Fig.  2  demonstrates  that  the  Raman  spectral  band  con¬ 
tours  as  computed  by  /)  quantum  mechanics,  using  Eq.  (Cl 4),  ii)  the  quantum  corrected 
correspondence  principle,  taking  the  classical  Ti— *0  limit  of  Eq.  (C14),  and  Hi )  the  quantum 
corrected  Newtonian  approach,  as  given  by  Eqs.  (3)  and  (5),  all  agree  closely  The  quantum 
spectrum,  suitably  broadened,  also  agrees  with  the  experimental  results  of  De  Santis  vt.  air" 
The  spectral  constants  for  the  quantum  and  correspondence  principle  calculations  are  from 
Huber  and  Herzberg-7  and  the  potential  energy  and  polarizability  are  as  described  in  Tables  I 
and  II  above.  In  this,  and  subsequent  figures,  the  integration  step  size  for  the  Newtonian 
molecular  dynamics  calculations  is  5x  10' w’ s. 

B.  Gas  phase  vibrational-rotational  spectra 

Figs.  3  and  4  show  the  fundamental  vibrational-rotational  isotropic  and  anisotropic  Raman 
bands  for  gas  phase  N,.  The  quantum  band  contours  are  calculated  from  Eqs.  (C15)-(C26).  the 
correspondence  principle  band  contours  are  from  the  classical  limit  of  Eqs  (CI5MC26),  and 
the  quantum  corrected  Newtonian  band  contours  are  from  Eqs  (4)  and  (5).  The  potential  and 
polarizability  functions  are  as  in  the  previous  section.  The  finite  time  step  used  to  obtain  the 
classical  trajectories  leads  to  a  small  offset  in  the  Newtonian  vibrational-rotational  calculations 
of  +5.5  cm  1  for  which  we  correct  before  presenting  the  spectra  The  ()  ( A 2 )  and  S 
(A7=+2)  branch  contours  are  seen  in  fig.  4  to  agree  well  among  the  quantum,  quantum 
corrected  correspondence  principle  classical,  and  quantum  corrected  Newtonian  approaches  We 
have  not  quantum  corrected  for  the  effect  of  anharmonicity  on  the  Q  (A./ •=())  line  shapes  and 
thus  while  the  Q  branch  areas  agree  among  the  three  calculation1,  the  line  shape  is  broader  in 
the  correspondence  and  Newtonian  classical  calculations  than  in  the  real  quantum  case,  as  is 
discussed  in  Appendix  E.  I  he  computed  quantum  Q  branch  agrees  with  ihe  experimental 
measurements  of  Bendtsen.17 

C  .  From  gas  to  liquid 

Having  demonstrated  that  the  quantum  corrected  classical  approach  can  provide  essentially 
correct  spectral  band  contours  for  gas-phase  rotational  and  vibrational-rotational  transitions, 
where  we  know  accurately  the  potential  energy  and  polarizability  functions  and  where  we  can 
compare  to  accurate  quantum  calculations,  we  now  turn  to  higher  densities,  where  out 
knowledge  of  the  potential  energy  and  polarizability  functions  is  less  certain,  and  our  ability  to 
compute  spectra  by  alternative  means  is  less  developed  We  will  treat  solutions  of  N  in  pro¬ 
gressively  higher  densities  of  Ar  at  298  K.  well  above  the  critical  temperature of  I5|  k.  and 
finally  in  liquid  Ar 

Work  of  others  which  relates  to  these  solution  studies  includes  ihat  of  Bratos  and 
Marechal''1  and  Bratos  and  l.eicknam.4'1  who  treated  Raman  band  shapes  in  solution  from  a 
correlation  function  viewpoint,  including  solvent  effects  in  a  stochastic  approach  Levesque. 
Weis,  and  ()xtobyJ!  have  treated  the  vibrational  dephasing  aspect  in  liquid  N\  and  knauss4-’  has 
studied  diatomics  in  solvent  Clarke  Miller,  and  Woodcock4,  have  computer!  liquid  N 
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HCi.  2  Rotational  anisotropic  Raman  spectra  lor  dilute  gas  phase  pure  N>  at  208  K  In  the 
upper  panel,  the  dashed  line  shows  the  gas  phase  computed  quantum  rotational  anisotropic 
spectrum  with  the  lines  broadened  so  that  they  can  be  compared  with  the  experimental  spec¬ 
trum.  shown  as  solid  circles,  of  De  Santis,  ci.  al.  in  Id  amagat  (20  atm)  N>.  The  experimental 
spectrum  is  given  in  arbitrary  cross  section  units  and  wc  have  normalized  ii  to  the  computed 
spectrum  The  solid  line  is  the  quantum  band  contour,  calculated  by  further  broadening  the 
lines  of  the  computed  spectrum  until  they  merge.  The  intensity  alternation  in  the  lines  is  Irom 
the  statistical  weights  of  the  nuclear  spin  states.  The  lower  panel  repeats  the  quantum  spectrum 
as  a  solid  line  and  compares  it  with  open  circles  showing  the  quantum  corrected  correspondence 
limit  classical  band  contour  and  open  triangles  showing  the  quantum  corrected  Newtonian  (clas¬ 
sical  molecular  dynamics,  classical  linear  response,  and  classical  ensemble  average)  spectrum 
The  Newtonian  spectrum  is  averaged  over  30000  single-molecule  24.6*10  L  s  time  histories, 
with  a  5*10  1,1  s  integration  step,  run  in  groups  of  300  non-interacting  molecules  The  isotro¬ 
pic  rotational  Raman  spectra  are  not  shown  as  they  are  zero 
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I  iCi  Fundamental  vibrational-rotational  Raman  spectra  for  dilute  gas  phase  N.  at  2')X  K 
Ihe  upper  panel  shows  the  isotropic  cross  section  and  the  lower  panel  shows  the  anisotropic 
cross  section  The  dashed  line  shows  the  computed  quantum  snectrum.  with  Ihe  lines  partialis 
broadened  Ihe  solid  line  is  the  quantum  band  contour  produced  by  further  broadening  ihe  in 
dividual  vibrational-rotational  lines  until  they  merge 
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MG.  4.  Fundamental  vibrational-rotational  Raman  spectral  band  contours  for  dilute  gas  phase 
N:  at  298  K.  The  solid  line  is  the  computed  quantum  band  contour,  from  Fig.  3.  The  open 
circles  show  the  quantum  corrected  correspondence  limit  classical  spectra  and  the  open  triangles 
show  the  quantum  corrected  Newtonian  spectra.  The  Newtonian  spectra  are  an  average  over  an 
ensemble  of  30000  single-molecule  24. 6x  10  v  s  time  histories 
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molecular  dynamics  for  rigid  molecules  and  compared  to  Raman  line  shapes.  Hanson  and 
McTague44-45  and  Frenkel  and  McTague46  have  computed  molecular  dynamics  for  N:  as  a  rigid 
dumbbell  and  have  used  this  approach  to  estimate  dipole  induced  dipole  effects. 

The  N3  intramolecular  potential  and  polarizability  tensor  are  as  before.  For  the  Ar-N, 
intermolecular  potential  we  use  a  sum  of  pairwise  atom-atom  potentials  between  the  Ar  and  N 
atoms,  each  of  the  form 

V(R)  =  A  exp[-2£(K-/?,.)]  -  2  tfexp|-£ (R-R,.)}  .  (10) 

in  which  ,4=5.793x10  J,  #=8.256x  10“22  J,  ^=!.51x]0111  m  and  R,  =2.94x10  m.  This 
potential,  which  would  be  a  Morse  potential  if  A  =  H,  is  due  to  Kistemaker  and  de  Vries,47  and 
has  been  shown  to  give  classical  trajectory  results47- 4X  in  reasonable  agreement  with  experimen¬ 
tal  rotational  relaxation  times47  for  N2  in  Ar.  The  Ar-Ar  interatomic  potential  is  the  Morse- 
spline-van  der  Waals  (MSV)  III  potential  of  Parson,  Siska  and  Lee.49  For  simplicity,  we  leave 
out  multi-body  effects  in  both  the  polarizability  function  (for  example  distortion  from  sphericity 
of  the  polarizability  of  the  Argon  atoms  during  collisions  and  distortion  of  the  N>  polarizability 
tensor)46-50  and  the  potential  function  (for  example  the  non-pairwise  additive  part  of  the  poten-  1 

tial)-M  which  should  be  included  in  a  more  refined  treatment.  The  Newtonian  approach  could 
also  be  used  to  study  these  multi-body  effects  (often  referred  to  as  collision-induced)  by  includ¬ 
ing  them  in  the  potential  and  polarizability  functions  and  comparing  the  results  with  experi¬ 
ment. 

To  facilitate  comparison  among  the  various  dilute  gas  to  dense  liquid  examples,  we  have 
not  made  the  appropriate  local  field  correction  which  depends  upon  the  refractive  index.57-  '1  To 
reduce  edge  effects  in  our  dynamics,  we  use  periodic  boundary  conditions,  replicating  a  trun¬ 
cated  octahedron  in  a  space-filling  solid  tessellation70  and  smoothly  reducing  (feathering)  our 
forces  to  zero  over  the  outer  0.01  nm  of  the  radius  of  the  inscribed  sphere  to  avoid  discontinui¬ 
ties  as  atoms  cross  the  boundaries.  Twenty  atoms  of  Ar  are  used  in  the  unit  cell.  The 
Newtonian  vibrational-rotational  spectra  are  again  corrected  for  the  small  +5  1  cm  1  shift 
caused  by  using  a  finite  5x10  16  s  time  step  in  our  trajectory  integration 

As  can  be  seen  in  Figs.  5  and  6,  there  is  drastic  change  is  spectral  line  shape  as  we 
proceed  from  dilute  gas  phase  to  a  density  of  1000  amagats.  (One  amagai  unit  of  density  is  the 
actual  concentration  of  the  particular  gas  at  0  C  and  1  atm  pressure,  which  would  be  2. 687x10-’' 
molecules  m  1  for  an  ideal  gas.)--14  The  spectral  line  shapes  for  rotation  (if  continued  through 
the  origin  to  also  show  the  Antislokes  scattering)  and  vibration-rotation  change  from  the  multi¬ 
ple  peaks  characteristic  of  gas  phase  free  rotation  to  the  single  peaks  characteristic  of  liquid 
phase  hindered  rotation  The  statistical  noise  caused  by  averaging  over  a  smaller  ensemble  of 
600  samples  can  be  seen  in  Fig.  5.  A  comparison  is  made  in  Fig.  6  with  the  experimental  O 
and  S  branch  band  shapes  measured  by  De  Santis  el.  al.*b  for  pure  N2  at  276  amagats,  and  the 
agreement  is  close,  even  though  the  solvent  is  Ni  instead  of  Ar,  indicating  that  the  line  shape  is 
not  particularly  sensitive  to  the  details  of  the  solute-solvent  intermolecular  potential  Note  the 
blue  shift  that  occurs  at  the  1000  amagat  density,  perhaps  explainable  by  an  effectively  steeper 
vibrational  potential  as  the  N2  is  crowded  up  against  the  repulsive  walls  of  the  solvent  The  evo¬ 
lution  of  these  spectral  band  shapes  with  density  can  be  compared  with  those  computed  for  the 
infrared  case  in  Paper  I,  and  with  the  example  Raman  band  shapes  computed  stochastically  by  ' 

Bratos  and  coworkers. 59  40 

D.  Liquid  phase 

We  compute  Newtonian  spectra  for  N2  dissolved  in  liquid  Ar  and  compare  with  the  exper¬ 
imental  anisotropic  spectra  of  Hanson  and  McTague.44-45 

In  Fig  7  we  show  the  rotational  spectra  The  Newtonian  calculations  are  made  at  the 
experimental  temperatures  with  the  experimental  densities  for  Ar  (792  amagats  at  83  K  and  749 
amagats  at  96  K),  with  N?  counted  as  another  Ar  atom  To  approximately  match  the  experi¬ 
mental  mole  fractions,  one  N,  is  run  with  28  Ar  atoms  in  the  unit  cell  at  83  K  and  one  N-.  with 
20  Ar  atoms  at  96  K,  although  no  direct  N2-N2  interaction  is  included  in  the  potential,  as  the 
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MG.  5.  Rotational  Raman  spectra  of  N>  in  different  densities  of  Ar  from  0  to  1000  amagats  at 
298  K  These  Newtonian  spectra  are  an  average  river  an  ensemble  of  MX)  time  histories,  each 
lasting  24.6k  10  |:  s,  for  a  system  consisting  of  a  single  N.  molecule  surrounded  by  20  Ar 
atoms,  using  periodic  boundary  conditions. 
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forces  are  feathered  to  zero  before  the  N?  in  the  adjoining  periodic  cell  would  have  been  felt 

The  agreement  of  quantum  corrected  Newtonian  calculation  and  experiment  is  seen  to  be 
quite  close,  but  not  exact.  It  is  not  known  how  the  blame  for  the  remaining  small  discrepancy 
should  be  apportioned  among  the  r )  inaccuracy  of  the  potential  and  polarizability  functions  we 
use,  in  particular  our  neglect  of  multi-body  (so  called  collision-induced  effects),4*’  51 1  u)  imper¬ 
fection  in  our  Newtonian  methodology.  Hi)  statistical  noise  due  to  inadequate  ensemble  sam¬ 
pling  with  only  200  runs,  and  n  )  possible  experimental  inaccuracy. 

I  or  the  vibrational-rotational  fundamental  band  shown  in  Fig  8,  the  agreement  with  ihe 
anisotropic  results  of  Hanson  and  McTague  is  excellent  except  near  the  very  peaks  of  the 
curves  This  disagreement  is  possibly  due  to  either  inaccuracy  in  our  theoretical  approach  or  to 
the  difficulty  which  Hanson44  cites  in  accurately  removing  the  much  stronger  isotropic  peak 
which  is  shown  in  the  upper  panel.  This  experimental  difficulty  is  compounded  by  the  factor  of 
one  tenth  in  Eq.  (7),  which  makes  the  ratio  of  parallel  to  perpendicular  experimental  intensities 
even  greater  than  the  ratio  of  isotropic  to  anisotropic  intensities  shown  in  Fig  8. 

IV.  CONCLUSION 

We  have  demonstrated  that  an  essentially  classical  mechanical  approach  can  be  used  to 
compute  the  Raman  spectra  of  dense  systems,  using  as  an  illustration  N-  in  compressed  Ar 
above  the  critical  point  (so  that  there  is  a  continuous  evolution  from  gas  to  liquid)  and  in  lower 
temperature  liquid  Ar.  Using  potential  energy  F(r)  and  polarizability  P(r)  functions  from  the 
literature,  with  no  adjustment  of  their  parameters,  we  have  achieved  close  agreement  with  the 
experimental  spectral  band  contours  from  dilute  gas  to  dense  liquid  In  addition,  we  find  very 
close  agreement  with  the  quantum  spectral  contours  for  the  dilute  gas  case  in  which  a  quantum 
solution  can  be  accurately  computed.  Since  the  band  shapes  are  derived  from  molecular 
dynamics,  we  can  have  confidence  that  the  computed  dynamics  provide  an  accurate  representa¬ 
tion.  within  the  limits  of  classical  mechanics,  of  the  actual  molecular  motions. 

As  with  the  molecular  dynamics  calculations  of  a  number  of  other  observables,  such  as 
transport  phenomena  and  chemical  reactions,  if  we  average  over  a  sufficient  number  of  quan¬ 
tum  states,  for  example  in  an  equilibrium  system  at  a  sufficiently  high  temperature,  so  that 
sharp  resonances  arc  not  observed,  then  measurables  are  usually  reasonably  well  reproduced  by 
classical  mechanical  calculations.  We  have  illustrated  two  related  ways  to  average  over  the 
resonances:  first,  lowering  the  resolution  of  observation  to  produce  the  band  contour  in  the  gas 
phase,  and  second,  broadening  the  resonances  by  perturbing  the  stationary  states  at  higher  den¬ 
sities  m  compressed  or  in  condensed  phases. 

The  results  in  this  paper  and  particularly  in  l  ig.  2  of  Paper  I  indicate  that  vibrational- 
rotational  band  shapes  can  be  quite  sensitive  to  the  precise  shape,  eg.  anharmonicity,  of  ihe 
intramolecular  potential,  but  that  neither  the  rotational  and  vibrational-rotational  band  shapes 
nor  the  corresponding  molecular  motions  arc  particularly  sensitive  to  (he  detailed  shape  of  the 
inlermolecular  potential  for  weakly  bound  solvent-solute  systems  Evidence  for  this  is  that  the 
band  shapes  are  well  reproduced  by  the  rather  crude  potential  forms  we  have  used,  and  that  the 
results  lor  Ar  versus  N\  as  a  solvent  do  not  appear  to  be  drasticallv  different,  as  shown  m  I  ig 
6 

t  hese  essentially  classical  calculations  agree  well  with  quantum  reality  for  translational  and 
rotational  motions  because  enough  states  are  accessible  lor  the  correspondence  principle  to 
apply  In  contrast,  for  N:  vibration  almost  the  whole  population  at  the  temperatures  considered 
is  in  the  ground  state,  and  thus  the  correspondence  principle  argument  is  not  valid  However, 
as  discussed  in  Paper  I.  molecular  vibrations  are  ordinarily  close  enough  to  harmonic  that  the 
special  parallels  between  classical  and  quantum  harmonic  oscillator  behavior''  'K  do  apply  and 
only  an  amplitude  correction  for  quantum  (here  zero-point)  versus  classical  thermal  motion  is 
needed.  The  anharmonicity  ol  the  oscillation  dries  not  enjoy  such  a  privileged  position,  and  we 
must  quantum  correct  our  classical  spectra  lor  two  anharnmnic  effects  if  we  wish  close  agree 
ment  with  experimental  realm  First,  a  small  shift  (the  order  of  a  percent  >  in  vibrational  line 
position  is  necessary  Second,  since  a  classical  anhaimonic  oscillator  vibrates  in  a  traiectory 
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MCi  M  Vibrational-rotational  Raman  spectra  of  N-  i n  liquid  At  The  Newtonian  spectra, 
shown  as  solid  lines,  are  computed  as  discussed  in  I  ig  7  1  he  tilled  circles  are  ihe  anisotropic 
experimental  results  of  Hanson  and  Me  I  ague  which  are  given  without  absolute  cross  section 
units  and  are  thus  normalized  to  the  theoretical  cross  sections 


which  is  no  longer  sinusoidal,  its  Fourier  transform  is  broadened  in  frequency  by  the  anhar- 
monicily.  In  contrast,  a  quantum  anharmonic  oscillator,  if  it  exhibits  essentially  only  a  single 
1—0  transition,  is  shifted  but  not  broadened.  Thus,  as  we've  seen,  in  the  very  narrow  Q 
branch  where  small  broadenings  are  observable,  an  additional  anharmonic  correction  would  be 
necessary  to  closely  match  the  quantum  width.  Where  anharmonicity  becomes  a  dominant  fac¬ 
tor,  an  improved  approximation  using  a  semiclassical  approach1  !  12  may  be  needed 

The  advantages  of  a  classical  approach  to  Raman  spectra  are  several:  /)  complex  molecular 
systems  can  be  treated  for  which  quantum  spectral  calculations  are  impractical,  for  example 
very  large  molecules,  clusters,  liquids,  solutions,  and  solids,  without  need  to  resort  to  normal 
mode  analysis,  //)  the  approach  can  be  extended  to  treat  time-dependent  and  non-equilibrium 
systems,  for  example  to  compute  the  Raman  spectra  during  the  course  of  chemical  reactions 
even  in  solution,  and  in)  experimental  spectra  can  be  analyzed  in  terms  of  atomic  motions  with 
the  help  of  classical  models  and  of  our  classical  intuition,  and  thus  understood  and  related  to 
molecular  dynamics  for  more  complex  molecules  and  systems  of  molecules.  In  fact,  the  quan¬ 
tum  correction  transformations  can  be  inverted  to  strip  away  the  quantum  aspects  from  an 
experimental  spectrum,  leaving  behind  a  classical  spectrum  which  is  truly  analyzable  from  a 
classical  view  point,  in  order  to  deduce  the  intra-  and  inter-molecular  motions  involved  in 
chemical  properties  and  processes.  One  further  step  toward  a  fully  classical  description  is  possi¬ 
ble,  as  Miller12  has  shown,  in  that  the  radiation  held  itself  can  be  treated  not  as  a  classical  field, 
as  we  have  implicitly  done,  but  rather  as  classical  harmonic  oscillators. 

In  conclusion,  the  greater  difficulties  in  computing  infrared  and  Raman  spectra  appear  to 
lie  not  on  the  side  of  nuclear  motion,  but  rather  on  the  electronic  side  of  the  Born- 
Oppenheimer  approximation,  i.e.  in  deriving  as  functions  of  nuclear  positions  sufficiently  accu¬ 
rately  the  potential  energy  L(r).  the  dipole  moment  #t(r),  and  the  polarizability  tensor  P(r) 
Given  these  functions,  classical  techniques  of  molecular  dynamics,  linear  response  and  ensem¬ 
ble  averaging  followed  by  simple  quantum  corrections  seem  capable  of  producing  quite  accurate 
infrared  and  Raman  spectra  for  dense  gases,  liquids  and  solids  (and  band  contours  for  dilute 
gases),  although  additional  quantum  corrections  may  well  to  be  needed  for  features  substan¬ 
tially  affected  by  vibrational  anharmonicity  such  as  overtone  and  combination  bands  and  other 
interactions  among  vibrational  modes  4  "4 
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APPENDIX  A:  LIGHT  SCATTERING  AND  POLARIZABILITY  TENSORS 

The  quantum  mechanical  expression  for  the  wavelength  weighted  Raman  differential 
scattering  cross  section  in  terms  of  energy  into  an  angular  frequency  range  </<»  and  a  solid 
angle  range  </U  is  given  in  the  polarizability  approximation  by1- 1 

X 4  -  5"  p  i<  /  i«  P-<"l/>i  (St<„,  (u )  (Alt 

(I  (oil  W 

in  which  27rA.  is  the  wavelength  of  the  scattered  radiation,  f>,  is  the  probability  that  the  system 
is  in  the  initial  state  /.  «  and  «  are  unit  vectors  in  the  directions  of  the  electric  vectors  ol  the 
incident  and  scattered  radiation,  respectively,  P  is  the  space  fixed  Cartesian  polarizability  tensor 
ol  the  system  ol  molecules,  and  to,  -  (A,—  I:  )/Tt,  where  /:.  and  /;  are  the  energies  of  the  final 
and  initial  states,  respectively 

The  (  artesian  polarizability  tensor  P  can  be  represented  as  the  .1x3  matrix  P  where  /  and 
/  are  each  one  of  the  space  fixed  axes  x,  v,  or  .-  This  tensor  can  be  decomposed  into  a  diago¬ 
nal  isotropic  tensor  and  a  traceless  anisotropic  tensor 
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p  =  p,„.  +  Pa,„s»  <A2) 

in  which 

P,w  =  {Tr  P  =  J  (A3) 


and 

P  =>P—  P  =F  —  a  I  (A4) 

1  utuso  1  r  iso  1  1  *  y.'i'vi 

I  being  the  unit  tensor  and  Tr  indicating  the  trace.  From  these  two  tensors  we  may  define  iwo 
polarizability  rotational  invariants,  and  which  measure  the  isotropy  and  anisotropy, 

respectively,  as 

}Tr(P„.P,J  =  «,L,  «  a2  (AS) 


and,  for  any  symmetric  polarizability  tensor, 

Tr[P /jniw P amsn )  Gannn  i  y  (A6) 

in  which  a  and  y  are  the  more  usual  spectroscopic  notation22  and 

o„„  =  +  py  +  p.)  ( A  7 ) 

. .  «  |((P„-Pl,)-,  +  (P,l-P:.-)‘>+(P,.--P«.)2  +  6(P>2+Pl’  +  P;-;)l  .  (A8I 


For  a  diatomic  molecule  with  polarizability  expressed  in  molecule  fixed  coordinates  with  the  : 
axis  along  the  internuclear  axis,  PU  =  P„.  and  all  the  off  diagonal  polarizability  components  are 
zero.  Thus,  for  a  diatomic,  Eqs.  (A7)  and  (A8)  in  molecule  fixed  coordinates  simplify  to 


=  j(P;  +  2P0)  ( A9) 

. -  |(Pc:  -  K)2-  ( A  1 0 ) 


We  can  write 


I*' 

k 


<l:tr 

dojciil 

d2tr 

daidil 


=  J,  f>,  I  </  |P,M,  l»  >  l:  hiw,,-U>) 

A<»  I.  f 

“  X,  P'  l</lPu».«.l,>l’  8(a*/,-w)  • 

am\o  /.  t 


(All) 

(AI2I 


the  isotropic  and  anisotropic  Raman  wavelength  weighted  cross  sections,  respectively 

In  order  to  avoid  the  awkward  properties  of  the  Cartesian  tensor  P  under  rotational 
transformations,  we  introduce  tensors  that  transform  properly/*"  We  construct  a  umt.ir> 
transformation,"1  represented  by  the  unitary  9x9  matrix  L.  the  components  of  which  are 
shown  in  Table  III 
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TABLE  111.  Components  of  the  U  matrix 
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U  transforms  the  Cartesian  tensor  P  when  written  as  the  9x  1  vector 

P  ~  \K  K  P  K  K  P,  h  K  K  I1  (Al 3) 

into  the  spherical  tensor  P.  T  indicates  the  transpose.  P  can  be  resolved  into  three  spherical 
irreducible  tensors 


p  _  p(OI  _)_  pill  +  pi  2) 

(A14) 

whose  components  P,'/'  are  for  P((l), 

P<T  =-4t(P.  +  Pi:i+P,)  . 

( A 1 5 ) 

for  Pn\ 

PtV  -  \lP^~K±VP:v~iPy:)) 

(Alb) 

Pi!"  =  ^Uh-iP„)  . 

<  A 1  7 ) 

and  for  P,2\ 

PW  «  !(PV,-P„  ±/Pv, ,±/P„) 

(A18) 

pi2,1  =  - ( ±  p,:  ±  p:v  +  /  p„  + 1  p,  i 

( A 1 9 ) 

Pi!:i  -  ^(2 P,.-Pv<-Pll) , 

(A20) 

as  can  be  verified  by  the  relations 

P  =  IP 

( A21 ) 

P  =  U  1  P 

(A22) 

in  which 

L'  1  =  (Ur>*  . 

(A2.ll 

the  asterisk  indicating  the  complex  conjugate  We  can  now  see  that  the  tensor  P  consists  of 
three  irreducible  tensorial  sets,  one  of  rank  0  with  I  component,  one  of  rank  1  with  3  com¬ 
ponents,  and  one  of  rank  2  with  5  components.  If  the  polarizability  tensor  for  the  system  is 
symmetric  (viz..  P„*»P„)  it  can  be  seen  from  Eqs  (A16)  and  (A17)  that  all  of  the  components 
of  P'",  the  irreducible  tensorial  set  of  rank  one,  are  zero. 

We  can  now  use  Eq.  (A22)  to  show  that 


-  12  - 


P,„  -  U"1  P(0’  -  U  1  0  0  0  0  0  0  0  0]1 

-  />r  I  =  {(K+K  +  K)  I 


( A24) 


in  agreement  with  Eq.  (A3)  and 


P 


anno 


=  D 1  P'2’  «  U  1  [0000  Pj-'  P,'2'  P,',-”  Pi] i*  P'>’]r 
=  P  -  U  1  P,(”  -  P  -  P„„ 


(A25) 


in  agreement  with  Eq.  (A4).  Using  Eqs.  (A24)  and  (A2S)  with  Eqs.  (All)  and  ( A 1 2 ) .  we 
have 


-7  Ip.  I</|/>00>  l/>l2  «(€u/f— oj)  (A26) 

/./ 

-Z>,  i  l</|P,i2,l'>l2S(^,-w)  .  (A27) 

anno  i.t  «i»-2 

These  two  equations  can  be  represented  by 

=  «,Ip,  i  \<i\PX'\f>\2  &(»»-*>)  (A28) 

i.t  /M—  / 

where  it  is  understood  that  y=  0  corresponds  to  the  isotropic  case  with  n,=  \  and  ./  =  2 
corresponds  to  the  anisotropic  case  with  1.  It  is  this  equation  that  will  provide  the  starling 
point  for  both  the  Newtonian  and  quantum  solutions  that  follow. 

In  addition,  it  is  often  convenient  to  expand  each  polarizability  componeni  in  a  Taylor's 
series  in  q ,  a  molecular  coordinate.  In  the  approximation  of  electrical  linearity,  we  keep  only 
the  first  two  terms,  yielding 


-Xs 


d2  c 


dojdil 
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dtodft 


4  d2cr 

'  du>d  n 


P  =  P"  4. 

rti  rn  ~ 


=  P"  + 

•n  ~ 
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<7 


(A29) 


Simularly,  this  can  be  extended  to  the  polarizability  invariants 

a,  =  «"  +  (A30) 

where  a,  is  for  j= 0  and  for  ./= 2. 


APPENDIX  B:  LINEAR  RESPONSE  THEORY 

The  linear  response  theory  of  Raman  spectroscopy  has  been  developed  by  Gordon.1' 14  1,1 
whose  approach  we  follow  here  and  also  discussed  by  others.17  -’9  h2  ^  The  first  step  is  to  con¬ 
vert  I.q  (A28)  from  the  Schrodinger  to  the  Heisenberg  picture  by  introducing  the  Fourier 
transform  of  the  Dirac  delta  function. 


ft(fti) 


2  n 
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e"“'  di  . 


From  Eqs.  (A28)  and  (Bl) 
d2 


-~7T  -  »,Ip,  t  \</W'\i>\2  y-  J  d,  exp 

du>di\  ,  '  ~  ,  2rr 


( E,  -  E, ) 


Because  |/>  and  |/>  are  eigenstates  of  the  system,  we  have 


1.17 


•ft  ft 


I  .  >ff  1  P  1 

l/>=e  "  !  i  > 


(BI ) 


/  i  .  (B2) 


(B3) 


and 
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,.i  it.i/1 1  ,.i  >H  i/ll 

<71  e  '  =</!<* 

in  which  H„  is  the  Hamiltonian  in  the  absence  of  the  radiation.  Therefore1'1 
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dmdil 


oo 

=  y ~f  dte-’-'n^P,  Z  <i\P,l‘)m[/><r\P:„"U)\i> 
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in  which 

P,:/’(r)  =  e'"“"n  P„V»( 0)  e  . 
By  closure, 

I  I/X/I  =  l  . 


(B4) 


(B5) 


(B6) 


(B7) 


and  therefore 

-4  d2a 


■x 


1  dudil 


oo 

=  2^  /  dt  £  p,  £  . 


Since  the  summation  overp,  is  just  an  ensemble  average  we  can  write 


16 


-X 


4  r/V 


27T 


oo 

$  dt  c  £  <P,!/'(0)P,!/'(r)> 


(B8) 


(B9) 


where  the  <>  now  indicates  an  ensemble  average.  Using  Eqs  (A5)-(A6)  and  ( A24 ) - ( A25 ) , 
Eq.  (B9)  can  be  transformed  into  Eqs.  (2)  and  (3)  of  the  text. 


APPENDIX  C:  QUANTUM  AND  CORRESPONDENCE  PRINCIPLE  SPECTRA 

Quantum  mechanical  expressions  for  the  Raman  cross  section  have  been  derived  using 
various  methods.22  61.63-66  We  will  follow  the  irreducible  tensorial  set60-62'67  approach  used  by 
Brodersen  and  coworkers68'7'1  which  is  applicable  to  many  molecular  geometries,  and  we  will 
use  a  somewhat  more  general  form  than  needed  for  our  'I  diatomic  in  an  attempt  to  clarify  the 
extension  to  more  complex  molecules.  Rigid  rotor  harmonic  oscillator  transition  probabilities 
will  be  used,  but  energies  (and  thus  frequencies  and  state  probabilities)  will  include  anhar- 
monic,  centrifugal  and  vibrational-rotational  coupling  terms. 

We  can  express  the  sum  of  the  tensorial  operator  matrix  elements  in  Eq.  ( A28 )  in  a  form 
which  applies  to  cither  the  isotropic  Q=0)  or  anisotropic  (J=2)  cases  as 

a  =  £  \<i\p,!;'U>r  .  (CD 

m  -  / 

If  we  now  assume  that  our  initial  and  final  states  are  separable  into  vibrational  and  rigid  sym¬ 
metric  top  rotational  states, 

|/>  =  \\JKM>  =  \\>\JKM>  (C  2) 

</ 1  -  CvV'AW/'l  =  <\’\<J’K'M'\  «  3) 

in  which  the  quantum  numbers  are  v  for  the  vibration,  J  for  the  total  angular  momentum.  K 
for  the  component  along  the  lop  axis,  and  M  for  the  component  along  the  space-fixed  r  axis, 
and  primes  indicate  the  final  state,  then  we  can  express  G  as 

G  «  £  \<v'\<J,K'M'\P^'\JKM>\v>\:  .  (C  4) 


If  the  coordinate  system  is  rotated  through  the  Euler  angles  >  .  V  to  coincide  with  the 
molecule-fixed  coordinates,  we  have69 

G  =  £  !  <  v'i  <  J  K  M  1  £  P/'1  D^'(-'P.-.x.-V)|yA,if>|v>|‘ 

*n  t  A  -  t 


( C  5 ) 
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in  which  D^,1  are  the  components  of  the  Wigner  rotation  matrix.71-72  The  components  of  A,''1 
are  the  same  as  those  given  by  Eqs.  ( A 1 5 ) - ( A20)  except  that  v,  y.  and  z  now  refer  to  the 
molecule-fixed  coordinates.  After  the  rotation,  ' 1  is  fixed  along  the  axis  of  vibration  and  can 
thus  be  extracted  from  the  rotational  matrix  element'’4 


G  =  £  I  £  <m'\P^'\m><J'K'M  ! /A', — X.  — ) |7AA/> I ’  .  (Co) 

ni-  /  k  -  ■  / 

The  absolute  value  squared  of  the  rotational  matrix  element  can  be  summed  over  M',  m  and  M 
symmetrically  using  Wigner  3j  symbols7-’ 

£  I  <J'K'MMK!n\JKM>  I2  =  (27+1X27+1  )(_^.  -(' 

M .  m.  \f 

-  (274  1)  C(7/7;AAA)2  (Cl i 

CUjJ  ;AAA )  is  the  Clebsch-Gordan  coefficient  for  the  coupling  of  angular  momentum  and  is 
tabulated  for  the  appropriate  cases  elsewhere.61’  It  vanishes  unless  the  selection  rules68 

A  -  A  =  A  A'  =  k  (C8) 

I j  -  /|  ^  r  ^  j  +  ./  tc4) 

and 

J'-J  even,  if  A'=A"=0  (CIO) 


are  fulfilled.  Thus  from  Fqs.  (A28),  (Cl),  (C6),  and  (C7), 


. % 4  d2<T 

'  doidil 


»,  £p,  *  |<v'|P'i'|v>|-'  (27+1)  CUjJ A AA )•’«(«.  -w)  . 


(Cl  I) 


in  which  a  degeneracy  factor  y  is  introduced  which  depends  on  the  symmetry  of  the  molecule 
Tor  our  'I  homonuclear  diatomic,  y  is  yJy  the  nuclear  spin  statistical  degeneracy.  Since  AA '•-=(). 
only  F.qs.  ( A 1 5 )  for  the  isotropic  and  ( A 20)  for  the  anisotropic  scattering  contribute  after  the 
polarizability  operators  have  been  rotated  to  internal  coordinates.  For  a  diatomic,  using  Fqs 
(A9)  and  (A10)  and  that  n.  is  '  for  /'= 0  and  I  for  /  =  2  we  can  replace  with  <*  and  1 

with  F.q.  (Cl  I )  now  becomes 


-r-*  d:<r 

dwdil 


£  p,  c /  <v"(i,  v>(  ( 2J  + 1 )  C  ( JiJ  .At))-  is (iu  - u) ) 


(CI2) 


with  the  isotropic  cross  section  being  given  by  /=  0  and  <»  ^<>  ...  and  the  anisotropic  by  /=-2  and 


«*  •=«».«, 

If  we  use  the  linear  expansion  of  «,  given  in  Fq  (A30U  the  vibrational  matrix  element 
for  harmonic  oscillator  wavefunctions  is  given  b>  4 


<v'|<»,"  +  a'  </  jv>  =  i*  ft,  ^  +  <>', 


THvfl) 

2aii  ,/M, 


ft.  I  V  • 


((13) 


in  which  a>„  is  the  angular  frequency  for  the  v+  I  — v  transition  and  m,  is  the  reduced  mass 

Using  Fqs  (C  l 2 1  and  (Cl 3)  and  the  Clebsch-Gordan  coefficient  C(J()J  .00)=  I  along  with 
the  values  of  the  other  Clebsch-Gordan  coefficients  given  in  fable  IV.  we  can  now  obtain 
expressions  for  the  the  pure  rotational  and  vibrational-rotational  differential  cross  sections  I  or 
a  pure  rotational  cross  section  we  have  Av-0  The  anisotropic  rotational  S  (A7^  +2)  branch  is 
given  by 


-.1  d  " 

'  do, dll 


£  p(v../i 


(ll  i im\it  ) 


3(7  +  1  U7  +  2)  .  AG. 7 +  2_+  /  (v.,/1 
2(27*31  *  ~n 


(a 


(CI4) 


The  isotropic  pure  rotational  cross  section  for  the  S  branch  is  zero 
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TABLE  IV.  Pertinent  squared  Clebsch-Gordan  coefficients 


A  J 

C(727';00)2 

-2 

37(7—1) 

2(27-1)  (27+1) 

0 

7(7+1) 

(27-0(27+3) 

+  2 

3(7+ 1 )  (7+2) 

2(27+0(27+3) 

For  a  fundamental  vibrational-rotational  transition  we  have  Av=I.  For  the  O  (A 7—2) 
branch,  the  isotropic  cross  section  is  zero  and  the  anisotropic  part  is  given  by 


d2  <T 
d(t)d{\ 

“  aniso 


£  p(v,7)  (a‘  amso 

v.J 


2  7T(v+D 

2<U()Wf 


37(7-1) 

2(27-1) 


£(v+17-2)-£(v,7) 

- — - CO 


(CIS) 


For  the  Q  (A7=0)  branch  we  have: 

4  B2!T 


I* 


dudil 

d2a 


dtodil 


=  X  P<v7)  («',«, )2 


=  £  p(v,y)  (a'ams0)2 

v.J 


7?(v+l) 

(27+1)  8 
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E(v+\,J)—E(v,J) 
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(06) 


7f(v+l) 


x  8 


2u)amr 

£(v+l,7)-£(v7> 


(27+07(7+0 

(27 —  1 )  (27+3) 


(Cl  7) 


And  for  the  S  (A7=+2)  branch 


*4  d2(T 
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2ooomf 

£(v+17+2)-£(v,7) 


3  (7+1)  (7+2) 
2(27+3) 
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(CI9) 


with  the  isotropic  cross  section  being  zero.  In  Eqs.  (C14MC18) 

-HlAv.J) 

p(v,y)  =  — —  - - 

t£gj  ^27 -F 1 ) 

v-o  y-o 

where  the  sums  in  Eqs.  (C14)-(C19)  are  over  all  the  initial  vibrational  and  rotational  states. 

The  final  approximation  is  to  use  energy  levels  evaluated  to  second  order  in  perturbation 
theory,  involving  the  third  and  fourth  derivatives  of  the  potential,  giving7* 

£(v7)  -  (j[o,(v+j)  +  £.7(7+1)  -  o,x,.(v+{)2  -  ae(v+\)JU+ 1)  -  D,7;(7+ l)2j  (C20) 
in  which 


B, 


If 


Ait  lr 

I,  =■  mrR,7 


(C21 ) 
(C22) 


In  the  above,  )"'(/?,.),  V"'{Rv)  and  V""(Rt.)  are  the  second,  third  and  fourth  derivatives  ol  ihe 
internuclear  potential  at  the  bond  length  R,  of  the  potential  minimum 

All  three  of  the  above  approximations  (linearity  of  the  polarizability,  rigid  rotor  harmonic 
oscillator  evaluation  of  the  transition  matrix-elements,  and  energy  level  evaluation  by  second 
order  perturbation  theory  through  third  and  fourth  derivatives  of  the  potential  function)  may  be 
extended  to  higher  terms  to  give  higher  accuracy.65-  7(>'7!i 

Using  these  formulas,  we  can  evaluate  the  quantum  spectrum,  and  then  obtain  the  band 
contours  by  broadening  the  8  functions,  for  example,  into  Gaussians 

8(<o/(— oj)  — *  (A/rrK'expl— A(<o„— <o)7]  ,  (C27) 

letting  A  decrease  until  the  individual  rotational  peaks  merge. 

By  evaluating  Eqs.  (C14MC26)  at  successively  lower  values  of  Planck's  constant,  until 
the  resulting  functions  converge  to  a  limit,  we  can  compute  the  correspondence  principle  classi¬ 
cal  differential  cross  section.  Our  Newtonian  calculations,  with  sufficient  runs  to  form  a  reliable 
ensemble  sample,  should  converge  to  these  correspondence  principle  results.  Slight  deviations 
can  result  from  the  approximations  in  Eqs.  (C14MC26).  but  they  should  be  negligible  in  the 
present  context. 


APPENDIX  D:  SEPARATION  OF  VIBRATION  AND  ROTATION 

We  have  already  made  the  approximation  of  separability  of  vibrational  and  rotational 
states  as  far  as  the  transition  matrix  elements.  If  we  make  the  further  approximation  of  separa¬ 
bility  for  the  energy  levels,  we  can  write  the  Raman  cross  sections  as  convolutions  of  vibra¬ 
tional  and  rotational  cross  sections.  If  we  use  Eq.  (Cll),  and  if  we  consider  a  vibrational- 
rotational  transition  corresponding  to  a  particular  A  A,  we  have 
Q 

-  n,  £,,,  I  <  v'  |  a*  |  v  >  | 2  ttj  (27+1)  CiJjJ'.K&K)2  .  (01) 

/  I 

in  which  the  superscript  Q  indicates  u  quantum  mechanical  differential  cross  section  and  /=() 

indicates  isotropic  and  7=2  anisotropic.  Let  l 

VIB(  i)  =  I  <  v'IPja  I v>  |7  (1)2) 

and 

ROT(  i)  =  n,  gj(2J+])CUjJ;K\K)2 .  (1)3) 

If  we  drop  any  coupling  terms  between  vibration  and  rotation  in  the  energy,  we  have 
Hwh  =  EW.J. K  )  -  E(x.J.K) 

=  [£(v')-E(v)J  +  [EU'.K')-EU.K)]  =  liw,  v  +  1iuJK  JK  (1)4)  j 

where  wv  v  and  wj  K  JK  are  the  angular  frequency  differences  between  the  vibrational  states  v’ 

and  v  and  the  rotational  states  7'A'  and  JK  respectively.  Then  we  can  separate  the  density  of  ! 


diudil 
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states  into  vibrational  and  rotational  factors 
p,  =  p(v,J.K)  =  p(v)pUA) 
in  which,  for  a  system  at  equilibrium. 


(05) 

(D6) 


pU.K)  = 


e-H  HJ.k) 


-/if  (7.  A)  ’ 


where  li=(ktiT)~\  kH  is  Boltzmann's  constant,  and  T  is  the  temperature.  The  angular  fre¬ 
quency  may  similarly  be  separated,  giving 

S (to /,  — to )  —  8(<ov  v+<o./a  to)  .  (1)8) 

Using  Eqs.  (D2)  and  (D3),  we  may  write  Eq.  (Dl)  as 

|x4  -  [I p(v)  u/fl(y)j  |£p<y.*)  /?ony)|  «<av  v+<*, ■*  Jh-U)  <iw 

This  can  be  expressed  as79 

(*4  ~~  -  [Ip(v)  P/fl(y>  8(«v-.v-w)|  *jIp(//0  KOTC./)  8(wyk  7A-to)| 


cy(y)  *  /?y0) 


in  which  the  wavelength  weighted  quantum  differential  cross  section  is  seen  to  be  a  convolution 
of  the  quantum  vibrational  cross  section  VQ(j)  and  the  quantum  rotational  cross  section. 
RU(.j).  Such  separability  approximations  are  common  in  correlation  function  analysis  of 
vibrational-rotational  spectra,*"  in  which  the  transformation  to  the  time  domain  can  be  made  by 
applying  the  frequency  convolution  theorem.27  (See  Paper  I  for  the  analogous  infrared  case.) 

APPENDIX  E:  QUANTUM  CORRECTIONS  TO  CLASSICAL  SPECTRA 

Our  original  wavelength  weighted  cross  section  given  in  Eq.  (Al),  and  our  separated  cross 
sections  W{j)  in  Eqs.  (D2)  and  (DIO)  and  RU(  j)  in  Eqs.  (C7),  (D3).  and  (DIO)  all  have 
the  basic  form 

.V^o)  =  S>,  |</|0|i>|2  »(«„-«)  (E-l) 

i.i 

in  which  S^fw)  is  the  quantum  spectrum  at  angular  frequency  tu,  /  is  the  initial  and  ./  the  final 
state,  p,  is  the  probability  of  state  /,  O  is  the  appropriate  operator,  and  w „-=(£, -£,)A.  in 
which  E,  is  the  final  and  E,  the  initial  state  energy.  If  we  substitute  -w  for  a>,  interchange  the 
subscripts  i  and /in  Eq.  (El),  and  recognize  that  |  </|Oji  >  |:=|  < / 10|/>  I2  and  fi(o>)=S(-i«>, 
we  then  find  that 

S^(—  w)  —  £  pi  |  <  / 10| / >  1 2  8(oi/,— o>)  .  (E2) 

i.i 

For  a  system  at  equilibrium,  p,=exp(j9ftai,>)p/,  and  thus 

“  exp {pfito,,)  -  exp(/37rw)  ,  (E3) 

which  is  just  an  expression  of  detailed  balance.  Because  of  the  6-functions  in  Eqs.  (El)  and 
(E2),  (u  can  replace  to,,.  The  correspondence  principle  classical  limit  of  Eq.  (E3)  is 
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-  lim  J^L  .  , .  (E4I 

Sc  (— oi)  n— ii  5^(— «u) 

in  which  the  superscript  C  indicates  classical.  Eqs.  (E3)  and  (E4)  can  be  understood  as  follows. 
Angular  frequency  to  connects  states  in  a  quantum  system  which  are  Hw  apart  and  of  appreci¬ 
ably  different  Boltzmann  probabilities  in  an  equilibrium  system,  and  thus  S^(cj)  and  Sy(- to) 
are  related  by  the  detailed  balance  factor  exp(/37?o>).  In  the  classical  correspondence  limit  as 
If— 0,  the  scattering  connects  states  of  only  infinitesimally  different  energy  7Fa>  and  thus  S(  { w) 
and  Sc  ( — to)  become  equal.  Comparison  of  Eqs.  (E3)  and  (E4)  suggest  that  one  should  con¬ 
sider  the  quantum  correction62-8182 

=  exp(/?W2)  (E5) 

•V  (to! 

which  symmetrically  and  simply  fulfills  the  requirements  of  both  equations.  The  factor  m/2 
inside  the  exponential  in  Eq.  (E5)  arises  because  S^f-w)  and  S0((o),  which  are  related  by 
detailed  balance,  lie  2to  apart  in  angular  frequency. 

Under  the  approximation  of  separation  of  vibration  and  rotation,  as  expressed  in  Eq 
(DIO),  we  will  quantum  correct  and  V^(j)  separately,  convoluting  the  results  where 

appropriate. 

1.  Rotational  correction.  From  Eq.  (E5)  above,  the  rotational  correction  is  just  the  detailed 
balance  factor,  and 

RU{j)  =  exp(/37/oj/2)  R(  (j)  .  (E6) 

Thus,  for  a  pure  rotational  band,  we  simply  multiply  the  classical  rotational  band  by  the  factor 
exp(07W2). 

2.  Vibrational  correction.  For  the  vibration  of  a  harmonic  oscillator  of  angular  frequency 
(on,  we  can  evaluate  the  quantum  vibrational  cross  section  V°(j)  in  Eq.  (DIO)  explicitly  as8* 

(*o>  -  £  (E7i 

v-0  V _ f  1/.Wr.  1  mr 


y°U)  - 1 


(P^)2  M<o„-co) 


■°  exp[~/3  ( )^a>0] 


2wmr  [  1  -exp(-/37?(o )  ] 


(PM)2  8((o„-co)  . 


V(  (  j),  the  correspondence  principle  7f— 0  classical  limit  of  Eq.  (E8).  is  given  by 

F<(./)  =  lim  VU(  j)  =  (2co‘V/3>  1  (Pi'A’)2  8((o,-(o)  (E9) 

1)  —41 

In  the  above  equations,  we  have  used  the  presence  of  the  8-function  to  exchange  w  for  a>(). 
From  the  ratio  of  Eqs.  (F.8)  and  (E9) 


Wj)  -  - ^T~=-T  Wj)  .  (E10) 

l-exp(-/37?w) 

giving  the  harmonic  oscillator  quantum  correction  for  vibration. 

For  a  vibrational-rotational  band,  by  Eq.  (DIO)  we  need  to  convolve  the  quantum  correc¬ 
tions  given  in  Eqs.  (E6)  and  (F.10).  The  result  is  that  to  quantum  correct  a  classically  calcu 
lated  vibrational-rotational  band  we  multiply  it  by 


1-exp  (-jSItu)) 

in  which  is  the  displacement  within  the  band  from  its  rotationless  center 

The  effect  of  potential  anharmonicity  is  essentially  the  same  as  presented  in  Paper  I.  The 
dominant  effect  is  to  cause  a  shift  in  the  classically  calculated  band  position  to  a  lower  quantum 
frequency  A  simple  correction  consists  of  evaluating  the  angular  frequency  for  the  rotationless 
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(v,7)  transition  (v+1,0)— (v.O)  first  quantum  mechanically  as  given  by  l:q.  (C  l  1 )  and  then 
classically  by  reducing  Planck's  constant  in  !:q  (('ll)  until  convergence  is  reached.  The  oil  set 
between  the  classical  and  quantum  rotationless  band  centers  can  then  be  used  to  shift  the  entire 
classical  vibrational-rotational  isotropic  and  anisotropic  band  contours  into  their  quantum 
corrected  position.  A  rougher  approximation.  -  4m  ,  (in  angular  frequency),  as  discussed  in 
Paper  I.  can  also  be  used. 

The  narrow  Q  branch  presents  a  special  problem  if  we  wish  to  examine  it  under  high  reso¬ 
lution.  as  there  is  an  additional  quantum  line  shape  correction  lor  anharmonicily  which  only 
becomes  important  for  it.  I  or  a  harmonic  oscillator,  the  quantum  and  classical  pure  vibrational 
spectra  are  the  same.  The  pure  vibrational  lines  of  a  quantum  anharntonic  oscillator  are  vari¬ 
ously  shifted  in  frequency  away  from  the  single  harmonic  value  by  the  anharmonicily.  but  not 
individually  broadened  In  contrast,  the  spectrum  of  a  classical  oscillator  is  broadened  by 
anharmonicily,  since  the  classical  oscillation  is  then  not  purely  sinusoidal,  and  in  addition  oscil¬ 
lators  with  different  energies  in  a  thermal  ensemble  have  different  frequency  distributions  If. 
as  in  the  case  treated  here,  Tia>„  for  the  vibrator  is  high  with  respect  to  A H  T so  that  only  a  single 
1  • — 0  quantum  transition  is  observed,  then  the  quantum  pure  vibrational  spectrum  will  be  a  sin¬ 
gle  narrow  line  shifted  in  frequency  by  the  anharmonicily,  while  the  classical  pure  vibrational 
spectrum  will  have  a  line  shape  considerably  broadened  by  anharmonicily  The  classical  line  will 
be  less  shifted  by  the  anharmonicily  because  the  classical  oscillator  samples  the  potential  at 
~-kH  /.  while  the  quantum  oscillator  samples  the  potential  at  the  considerably  higher  energy 
— 77 ojt i  where  the  effect  of  anharmonicily  is  greater.  The  anharmonic  effect  on  the  width  of  the 
classical  and  quantum  pure  vibrational  spectra  is  ordinarily  (/.<•.  in  ().  P,  R.  S  branches) 
overwhelmed  by  the  much  larger  width  produced  by  rotational  broadening  for  the  0  branch, 
the  rotational  broadening  collapses  to  a  narrow  width,  the  anharmonic  broadening  of  the  classi¬ 
cal  vibrational  line  shape  becomes  apparent,  and  a  quantum  line-narrowing  correction  for  v  ibra¬ 
tional  line  width  must  be  made  if  the  quantum  Q  branch  line  shape,  and  not  just  area,  is  to  be 
reproduced  by  a  classical  calculation.  This  could  perhaps  be  accomplished  by  deconvoluting  the 
vibrational  width  from  the  classical  Q  branch  The  true  quantum  width  of  the  1—0  Q  branch 
can  also  be  more  closely  approximated  classically  by  the  artifice  of  computing  the  spectrum  for 
a  classical  harmonic  oscillator,  in  which  the  pure  vibrational  line  width  (ignoring  collisional  and 
radiation  damping  effects)  collapses  to  a  ft-funeiion.  and  the  remaining  width  is  from 
rotational-vibrational  coupling,  and  then  shifting  the  line  in  angular  frequency  by  -  4nr  v.  to 
correct  for  the  anharmonic  shift  The  detailed  rotational  shape  of  the  Q  branch,  however, 
would  not  be  correct,  as  the  rotational-vibrational  coupling,  quantum  mechanically  the  <«,  in 
iiq.  IC25>,  will  be  changed  by  omitting  the  anharmonicily  The  sensitivity  of  the  shape  of 
actual  (as  contrasted  to  theoretical  rigid-rotor  harmonic  oscillator)  Q  branches  to  anharmonicily 
through  vibrational-rotational  coupling  has  often  been  insufficiently  emphasized  in  attempts  to 
assign  only  a  vibrational  origin  to  isotropic  line  shapes  In  conclusion,  ii  may  be  asking  too 
much  of  classical  mechanics  to  exactly  reproduce  the  shape  of  a  very  narrow  quantum  line  like 
the  Q  branch  in  a  system  in  which  there  are  the  twin  obstacles  to  classical  convergence  of  i) 
population  in  only  a  single  vibrational  quantum  state  and  u)  a  major  role  played  by  anharmom- 
city  Perhaps  the  introduction  of  aspects  of  quantum  mechanics  directly  into  the  molecular 
dynamics,  in  other  words  a  senuclaxsical  approach.1  1  is  the  better  solution  for  this  case  In 
the  correspondence  principle  and  Newtonian  classical  spectra  presented  here  we  have  not  made 
any  anharmonicily  quantum  correction  to  the  Q  branch  line  shapes  and  thus,  while  the  0 
branch  areas  are  essentially  correct,  the  shapes  of  the  correspondence  principle  and  Newtonian 
0  branches  are  somewhat  distorted  from  quantum  reality,  as  can  be  seen  m  lag  4 
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